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Abstract. We consider algorithms for finding the optimal location of a simple transportation device, 
that we call a moving walkway, consisting of a pair of points in the plane between which the travel speed 
is high. More specifically, one can travel from one endpoint of the walkway to the other at speed v > 1, 
but can only travel at unit speed between any other pair of points. The travel time between any two points 
in the plane is the minimum between the actual geometric distance, and the time needed to go from one 
point to the other using the walkway. A location for a walkway is said to be optimal with respect to a 
given set of points if it minimizes the maximum travel time between any two points of the set. We give a 
simple linear-time algorithm for finding an optimal location in the case where the points are on a line. We 
also give an Q(n log n) lower bound for the problem of computing the travel time diameter of a set of 
points on a line with respect to a given walkway. Then we describe a 0(n log n) algorithm for locating 
a walkway with the additional restriction that the walkway must be horizontal. This algorithm is based 
on a recent generic method for solving quasiconvex programs with implicitly defined constraints. It is 
used in a (1 + e) -approximation algorithm for optimal location of a walkway of arbitrary orientation. 
Finally, we consider several variants of the transportation model, including escalators, used for traveling 
between two floors of a building, and elevators, which are strictly vertical. 



1 Introduction 

Finding the optimal location of a facility modeled as a geometric object is an old, well-studied topic in 
computational geometry. One of the simplest instances of such a problem is finding the smallest disk cov- 
ering a given set of points, which can be done in linear time, or its min-sum counterpart, the Fermat- Weber 
problem, which is difficult to solve exactly in usual models of computation. Surveys on both combinatorial 
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and geometric facility location problems can be found in [1 1, 12, 20]. The field is also tightly related to the 
theoretical notions of center, median and depth of a set of points [6]. 

The optimality of a location is often defined with respect to a distance measure. Recently, an interest 
has emerged in the definition and analysis of time distances, measuring the travel time between points in 
the plane given a geometric transportation facility, such as a highway or a geometric network made of line 
segments [1-3, 5, 16]. In general, the transportation facility is a subset of the plane where the travel speed is 
high, and the travel time is naturally defined as the minimum time needed to reach a point. Primary concerns 
in previous works on these notions include in particular efficient constructions of the Voronoi diagrams for 
these distances. 

In this work, we consider optimal location of geometric transportation facilities on a line and in the 
plane. The transportation facilities we define are probably the simplest one can think about, and consist of 
a pair of points between which one can travel faster than anywhere else in the plane. This simple definition 
can be used to model for instance moving walkways in airports, shuttle connections or tunnels in a city. For 
simplicity, we refer to these as moving walkways. Given a set of points and a moving walkway, we define 
the travel time diameter as the maximum travel time between two points in the set. A moving walkway is 
said to be optimal if it minimizes the travel time diameter. 

Moving walkways on a line are defined and analyzed in section 2. We first provide an fl{n log n) lower 
bound on the computation of the travel-time diameter for a given walkway in the algebraic computation tree 
model. Then we give a simple linear-time algorithm for finding the optimal moving walkway on a line. 

Section 3 summarizes some previously known results on LP -type and quasiconvex programs. 

In section 4, we define moving walkways in the plane, and provide an 0(n log n) deterministic algorithm 
for the diameter. We show, using Chan's technique for implicit quasiconvex programs, that this implies an 
0(n log n) algorithm for finding an optimal horizontal walkway. We then derive an algorithm for computing 
a (1 + ^-approximation of an optimal moving walkway with arbitrary orientation and constant speed in time 

0(e~ 1 nlogn). 

Section 5 contains a discussion on how the transportation model can be restricted to yield the escalator 
and the elevator models, and on the complexity of the corresponding location algorithms. We also give an 
algorithm for computing the travel time diameter of a set of red and blue points with respect to k elevators. 

Related works. We are not aware of many previous contributions on the topic of optimal transportation 
facility location based on time distances. 

In [8], we presented a family of geometric optimization problems in which the objective function was in 
a min-max-min form. As an example, we proposed the problem of finding the location of a moving walkway 
on a line such that the maximum travel time between n given pairs of source-destination points is minimized. 
The current paper presents extensions and generalizations of this result. 

Recently, an effort has been made by several authors to solve the problem of locating a highway that 
minimizes the maximum travel time among a set of points [4]. In this work, a highway is modeled as a line 
in the plane on which the travel speed is higher, and customers can enter and exit the highway at any point on 
the line. In particular, they give an algorithm for the case where the travel speed on the highway is infinite, 
and another algorithm for the case where the highway is restricted to be vertical. A 0(n 2 logn) algorithm 
for the general problem is also given. Various other models are studied by Korman and Tokuyama [18]. The 
goal of this work is similar, but the transportation model that we study is different. 



2 Moving walkways on a line 



Suppose a new moving walkway is to be installed in a long corridor, for instance in the concourse of an 
airport. The moving walkway is used by customers to go from source to destination points in the corridor. We 
model the corridor as the real line, source and destination points as real numbers, and the moving walkway 
by an interval [a, b] on the line. We denote by v > 1 the speed on the moving walkway, and assume that the 
speed outside the moving walkway is 1 . Customers can only enter the walkway at point a, step down of it at 
point b, and always follow the shortest path. This is illustrated in figure 1 . 
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Fig. 1. Moving walkways on a line. The travel time between two points is the minimum between a) the direct distance, 
b) the travel time using the walkway. 

The travel time between a point s of the line and any other point t > s, given a walkway [a, 6], is a time 
distance [2, 5], defined as follows: 

rftimc(s, t, a, b) = min{i — s, Is — a\ + \t — b\ H — (b — a)}. 

v 

Note that by assuming t > s we implicitly make the moving walkway bidirectional: for going from t to 
s, we enter the walkway at point b and exit at point a, which makes the travel time a symmetric distance 
function. 

We consider the computation of the travel-time diameter and the optimal location of such a moving 
walkway. 

2.1 Travel time diameter 

A first problem is, given a set of points and a walkway [a, b], to compute the travel time diameter, defined as 
the maximum travel time between two points of the set. 

Problem 1 Given a set P of n real numbers, a real number v > 1, and a pair (a, b) of real numbers, find 

max dt- ime (s,t,a, b). 

s,t£P:s<t 

Theorem 1. The travel time diameter for a moving walkway on a line can be computed in 0(n log n) time. 

Proof. Without loss of generality, we assume that max{i e P} = 1 and min{a; € P} = 0. We distinguish 
the four intervals [0, a], (a, (a + b)/2], ((a + b) /2, b] and (b, 1], and call these I, II, III and IV, respectively. 



We remark that if a pair s, t is not separated by the point (a + b)/2, then the travel time between s and t 
is \t — s\, i.e. we never use the walkway. From this, the pairs of points between which the moving walkway 
is used for traveling can only be of four types: II-III, I-III, II-IV and I-IV. For all the other cases, the travel 
time is equal to the maximum absolute difference within the corresponding sets of points. 



Fig. 2. Computing the travel time diameter for points in the interval [a, b] amounts to computing the angular diameter of 
points on a circle (here v = +00). 

For the case II-III the computation of the diameter is essentially equivalent to computing the maximum 
angle between two points on a circle (see figure 2). This can be achieved in linear time once the input is 
sorted using rotating calipers [23]. 

The case I-III is solved by computing the farthest point from a on the same circle, which can also be 
achieved by scanning the sorted array. The case II-IV is symmetric to the case I-III. Finally we know that 
all pairs of type I-IV use the walkway, so the largest travel time is dtime(0, 1, a, b), which can be found in 
constant time. □ 

Theorem 2. There is an J7(nlogn) lower bound on the computation of the travel time diameter for a 
moving walkway on a line in the algebraic computation tree model, even when v = +00. 

Proof. The proof of the lower bound is by reduction from the problem of set disjointness, consisting of 
checking whether A n B — for two sets of numbers A and B. The proof is similar to the proof of the same 
lower bound for the problem of computing the Euclidean diameter of a set of points. 

We let v = +00, a = 0, b = 2 and scale the elements of A and B so that they are all strictly smaller than 
1. We then let 



Then the travel diameter of the points in P with respect to the moving walkway a, 6 is equal to 1 iff A n B ^ 
0. Note that since 1 is an upper bound on the travel time diameter, the lower bound also holds for the decision 



2.2 Minimum travel time diameter location 

We now consider the following moving walkway location problem. 

Problem 2 Given a set P of n real numbers, a real number v > 1, find a moving walkway [a, b] of speed v 
that minimizes the travel time diameter of P: 



a 



b 




P = AU {x + 1 : x e B}. 



version of the problem. 



□ 



min max dtime(s, t, a, b). 

a,b s,t<£P:s<t 



There is a very simple deterministic algorithm solving this problem in 0(n) time. It is a surprising 
observation that we are able to find a moving walkway that minimizes the diameter faster than we are able 
to compute the diameter for a given walkway. This observation has already been made for the construction 
of optimal vertical highways [4]. We assume w.l.o.g. that minja; G P} = and max{i G P} = 1, and 
denote by OPT the optimum value min a ^ max s j 6 p :s < t d t i mc (s, t, a, b). 

Lemma 1. 

v 

OPT < 



2v- 1 

Proof. We define a moving walkway by a = (1 - ^=i)/2 = ^ and b = + l)/2 = f^. Let 

ci = 1 — 2^rj be the symmetric of point with respect to a. Similarly let c 2 = ^xhi be the symmetric of 
1 with respect to b. Then going from any point in [0, ci] to a point in [c 2 , 1] using the walkway cannot take 
more time than going from to 1 . 

To go from to 1, the time needed is 

o+(l-6)+ -(b-a) = '" 



v v ' 2v - 1 

For every pair in [0, c 2 ], the time distance is at most (the length of the interval). The same observation 
holds for every pair in [ci , 1] . □ 

In what follows, we assume we are given a set of points P and we try to find properties of its optimal 
moving walkway [a, b]. Let s be the smallest point in P such that, to go from to s, it is better to use the 
moving walkway. We know that s < b, as to go from to any point larger than or equal to b it is always 
better to use the moving walkway. Symmetrically, let r be the largest point in P such that to go from 1 to r, 
it is better to use the moving walkway; we know that r > a. 

The following facts are easy to verify: 

Fact 1 For every pair (p, q) enclosing the walkway, i.e., p < a < b < q, the walkway is used. 

Fact 2 If we have p < p' < a, and the walkway is used for the pair (p' , q), it is also used for the pair (p, q). 

Fact 3 If the walkway is used for the pair (p, q) (with p < q), and p' is between p and a, and q' is between 
q and b, then the walkway is used for the pair (p' , q'). 

Lemma 2. Let P be a set of points and [a, b] an optimal moving walkway. The moving walkway is used only 
for pairs (p, q) where p G [0, r] and q G [s, 1]. For every other pair, the facility is not needed. 

Proof. We first show that for a pair (p, q) in the interval [0, s), the facility is never used. We distinguish three 
cases. 

1. if (p, q) belongs to [0, a], the walkway is trivially not used. 

2. if p < a and q belongs to [a, s), then using the walkway for going from p to q implies that we should 
use it for going from to q (see Fact 2). This contradicts the definition of s. 

3. if both p and q belong to [a, s), then using the walkway for going from p to q implies that we should use 
it also for going from p' to q, where p' is the symmetric of p with respect to a. We are then back to case 
2, leading to a contradiction. 



In conclusion, the walkway is never used for pairs in the interval [0, s). The same reasoning holds for pairs 
in(r, 1]. □ 



Lemma 3. The optimal moving walkway is such that a = | and b = ^il. 

Proof. We have not proved that all pairs (p, q) with p£ [0, r] and q e [s, 1] require the use of the moving 
walkway. What we know however is that the moving walkway will be used for the pairs (0, 1), (0, s) and 
(r, 1). What we aim to show is that finding the optimal moving walkway is equivalent to finding the optimal 
moving walkway for these pairs. 

First note that the optimal walkway for (0, 1), (0, s) and (r, 1) is such that a = | and b = ^p. More- 
over, it implies that d t i me (0, 1, a, b) = d t i me (0, s, a, b) = d timo (r, 1, a, b). We also have d t i me (0, 1, a, 6) > 
dtime(^) s 7 &)> meaning that we also optimize the distance for (r, s). 

For every pair with p € [0, r] and q € [s, 1], we have one of the following: 

- if p < a and q > b, the walkway is always used (Fact 1), 

- if p < a and q < b, the walkway is always used since it is used for (0, s) (Fact 3), 

- symmetrically, if p > a and q > b, the walkway is used since it is used for (r, 1), 

- if p > a and q < b and the walkway is not used, then the time distance for (p, q) is no greater than the 
shortest distance from p to q using the walkway, which is no greater than the time distance for (r, s) . If 
it is used, it is no greater than the time distance for (r, s) . 

So, in every case, the time distance is upper bounded by that of (0, 1), (0, s), (r, 1) or (r, s), the maximum 
of which is minimized with the solution described above. □ 

This lemma showed that the optimal moving walkway is uniquely determined by the points r and s. 

Lemma 4. Let 

- n be the largest point of P whose value is less than or equal to 1 — T^ry 

- s\ be the smallest point of P whose value is more than or equal to ZiilLzll+Jl+l 

- T2 be the largest point of P whose value is less than or equal to ^~y~g~^ 

- S2 be the smallest point of P whose value is more than or equal to 2 v-i - 

Then, either r = r\ and s = Si or r = r^. and s = S2. 

Proof. By Lemma 1 and Lemma 3 we know that either r < 1 — 2 v-i or s — i -l ' ^ or otherwise the 
walkway (a = r/2, b = leads to a diameter greater than ^xh\ • ^ e treat tne ^ rst case ' t ' le otrier one 
being symmetric. 

If r < 1 — 2v-i ' we P rove tn at r = r\ and s = s\. We proceed in two steps: first we will show that 
r = n is optimal, and then prove that s = s\ independently of the choice of r. 

By contradiction, imagine that the optimal r is not r\. Because we assumed that r < 1 — > this 
means that there exists a point q between r and 1 — 2 v _ 1 , as r is not the largest. Since q > r, and r is, by 
definition, the largest point for which we should use the walkway to go to 1, we know that the walkway is 
not used from q to 1, leading to a diameter greater than 2^3T> an ^ contradicting Lemma 1. To conclude, we 
know that in the case where r < 1 — 2 v _ 1 , r is the largest point of P satisfying that constraint. 

Now we prove that s = si. By definition of s, it is the first point in P such that, to go from to s, it is 
better to use the walkway. Thus we obtain the following equation: 




s > 



r(v — 1) + v + 1 



3v-l 



By definition, s is the smallest point in P satisfying that constraint, and thus, as claimed in the statement 
of the Lemma, s = s\. □ 



Theorem 3. Problem 2 can be solved in 0(ri) deterministic time. 

Proof. By Lemma 4, we know that there are two possibilities. So first, we determine the values of n, s\ 
and s 2 in linear time (we scan the set twice). By Lemma 3, the optimal moving walkway is given by a = ~ 
and b = ^±1, so we have two candidates. 

By Lemma 2 and the arguments developed in the proof of Lemma 3, we know that the largest time 
distance for pairs (p, q) with p E [0, r] and q £ [s, 1] is achieved by one of the pairs (0, 1), (0, s), (r, 1) 
or (r, s). The largest time distance for pairs where both points lie in [r, s] only depends on the smallest and 
largest point in that ranges and can also be found by scanning the set. For the pairs in [0, r] or in [s, 1], the 
worst distance achieved is max (r, 1 — s). 

In other words, we compute the worst time distance for both candidates (for pairs using and not using 
the walkway), and keep the best. □ 



1/2 s 



Fig. 3. Illustration of the moving walkway location algorithm on a line. 



3 Randomized optimization algorithms 

In this section, we briefly recall useful results on generic optimization problems. These results will be used 
in the next section on moving walkway location in the plane. 

The problems we consider are constrained minimization problem of the form 

min^ y 
subject to y > fi(x) Vi = 1,2, ... ,n, 

where x G y € R and {/.;} is a collection of continuous functions of the form /, : W.^ 1 BL In 

order to have a unique minimum, we can define the min operator as a lexicographic minimum on the vector 
made of the value y and the coordinates of x. Geometrically, the problem can be interpreted as seeking a 
lowest point on the upper envelope of a set of (d — 1) -dimensional surfaces fi(x). If these surfaces are 
hyperplanes, then we have a linear programming problem. If they are defined by convex functions, then we 
have a convex programming problem. In both cases, there exist efficient randomized algorithms running in 
0(n) expected time. 

In fact these algorithms do not apply solely to linear and convex programming problems, but to a well- 
defined family of so-called LP-type problems [22]. An abstract c?-dimensional LP-type problem is a collec- 
tion of constraints H and a function w : 2 n i— > W associating a value in a totally ordered set W to each 
subset of H. The goal is to compute w(TL). Moreover, w has to satisfy the following, for any H C H and 

heH: 



1. w(H) — w(B) for some B C H of size at most d; B is called a for H, 

2. u;(i?U{fe}) > tu(ff), 

3. If w{H) = w{B) and B C H, then U {h}) = w(H) ^> w(B U {h}) = w(B). 

Provided that d is a constant and the computation model allows the execution of some elementary operations 
in constant time, then an abstract LP-type problem can be solved in O(n) randomized expected time. 
In the above problem, the set W is the set K, a constraint h is a surface fi(x), and the value w(H) 
of a set H of constraint is the height of the lowest point above all the surfaces in H. If the functions fa 
are convex, the dimension of the LP-type problem is d, since an optimum w(H) is defined by d surfaces in H. 

A function / is called quasiconvex if all its (lower) levels are convex, i.e. if the set {x : f(x) < t} is con- 
vex for all value of t. Alternatively, a quasiconvex function is such that f((x + x') /2) < max{f(x),f(x')} 
for any pair x, x'. Quasiconvex programming is essentially the problem defined above in the case where all 
the functions /, are quasiconvex [14, 10]. Amenta, Bern and Eppstein showed the following: 

Theorem 4 ([7]). Any quasiconvex program of dimension d forms an LP-type problem of dimension at most 
2d+l. 

As a direct consequence, any quasiconvex program with n constraints can be solved in O(n) randomized 
expected time. 

Recently, Chan introduced a generic technique for solving some LP-type problems with a large number 
of constraints [9]. He motivated the introduction of this method by the problem of computing the Tukey 
depth of a set of points. It has since been used in a number of other contributions, including a new algorithm 
from Morin for computing the Zonoid depth of a set of points [19], and a facility location algorithm from 
Eppstein and Wortman [15]. We give the main result applied to quasiconvex programming, using notations 
from [14]. 

Theorem 5 ([9]). Let V be a set of input values and H a set of quasiconvex functions. If there exists a 
function g : 2 V i— > 2 H such that: 

1. the quasiconvex program defined by the constraints g(B) with B C V can be solved in O(l) time 
whenever \B\ = 0(1), 

2. there exists an algorithm running in 0(D(\P\)) time for deciding whether a given pair (x,y) is a 
feasible solution for the program defined by P C V , i.e. whether y > f(x) V/ € g{P), and such that 
D(n)/n £ is monotone increasing for a constant e > 0, 

3. there exist two constants a < 1 andr G N such that we can, in time 0(D(\P\)), decompose the problem 
P C V into r subproblems P\ , Pi, . . . , P r of size at most a\P\ and such that g(P) — U[=i 9 (Pi) 

then for any P € V , we can solve the quasiconvex program defined by g(P) in randomized expected time 
0(D{\P\)). 

So under some simple conditions, we are able to solve a quasiconvex program in essentially the same time 
as needed to decide whether a given solution (x, y) is feasible. This is achieved by a generalization of the 
standard randomized incremental algorithm in which the bases are evaluated recursively, and the satisfaction- 
violation test is handled by the algorithm provided by requirement 2. 

This is interesting in situations where the input is not an explicit list of constraints, but rather some data 
P that can be mapped to a set of constraints g(P) that is much larger. This yields what are called implicit 
quasiconvex programs. This is particularly useful for our minimum-diameter problems, since the diameter 



is defined by a pair of points in a n-point set, yielding (™) constraints. But if the above conditions are met 
and we can test the diameter for a particular solution in time D(n), then we can find the minimum-diameter 
solution in time 0(D(n)) as well. 

4 Moving walkways in the plane 

In this section, we generalize the previous moving walkway model to the plane. A moving walkway is 
defined by two points a and b, and customers can go from a point s to a point t either by going from s 
to t directly, or by entering the walkway at one of its endpoint and exiting at the other. The walkway is 
parameterized by a real number v so that the time to go from a to b is -d(a, b). The distance measure d(., .) 
used in the definition of the travel time must be convex and satisfy the triangular inequality. For simplicity 
we assume it is the Euclidean distance, keeping in mind that all the following can be generalized to other 
distance measures. Note that the moving walkways are bidirectional, hence the travel time is a symmetric 
distance function: 

^timc( s 7 1, a, b) = min{rf(s, t),d(s, a) + -d(a, b) + d(b, t),d(s, b) + -d(a, b) + d(t, a)}. 

v v 

We first solve a location problem with n source-destination pairs, under the restriction that the moving 
walkway must be horizontal. Then we give an algorithm for computing the travel time diameter for a given 
moving walkway (not necessarily horizontal), and from this derive an algorithm for finding optimal hori- 
zontal walkways. We show that this algorithm can be used for obtaining an approximation algorithm for 
locating moving walkways with arbitrary orientation. 

4.1 Horizontal moving walkway location for n source-destination pairs 

As a preliminary, we consider the problem of minimizing the maximum travel time between a number of 
designated source-destination pairs, with the restriction that the moving walkway must be horizontal. We let 
(Px,P y ) denote the x and y-coordinates of a point p. 

Problem 3 Given n pairs {(si, fj)}™ =1 of points, and a real number v > 1, find a horizontal moving walk- 
way [a* , b*) of speed v minimizing the maximum travel time between a pair: 

min maxd timc (s l ,U,a,b), 

a,b:a y —b y i 

The following lemma states that if one travels from left to right, then the moving walkway, if it is used, 
will be used from left to right. 

Lemma 5. Given a pair of points a, b such that a y = b y and a x < b x , and a pair of points s, t such that 
s x < tx, we have: 

rftimc(s,t,a, b) = min{d(s,t),d(s, a) + -d(a,b) + d(b,t)}. 

Proof. As it is the case on the line, pairs of points s, t for which the walkway is useful are separated by the 
bisector of the segment [a, b]. Then we can see that if we use the walkway from b to a, using it from a to b 
can only reduce the travel time. □ 



Lemma 6. Problem 3 is LP-type. 



Proof. We prove that the functions 

fi(a, b) = d t i mc (si,ti,a,b) 

are quasiconvex. From Theorem 4, this implies that the problem is LP-type. The domain of these functions is 
M 4 , however we only consider their restriction to the subspace of points a, b such that a x < b x and a y = b y , 
for which quasiconvexity is preserved. We also assume, without loss of generality, that s ix < t ix . 
We define the following function g. L , defined on the same domain as ff. 

gAa, b) = d(si, a) H — d(a, b) + d(b, tA. 
v 

This function is convex, since it is a positively weighted sum of convex functions. 

Now we have /j(a, b) — min{d(si, it), <ft(a, b)}. The value d(sj, ti) does not depend on a and 6, and 
can therefore be considered as a constant. So /; is defined as the minimum between a constant and a convex 
function. All the levels of such a function are necessarily convex (see figure 4). □ 

LP-type problems are known to be solvable in expected time linear in the number of constraints [22], 
which implies that Problem 3 is solvable in 0(n) expected time. 




Fig. 4. Functions defined as the minimum between a convex function and a constant are quasiconvex. 



Remark 1. The restriction on the direction of the moving walkway is useful in making the functions ft 
quasiconvex. If this assumption is lifted, then it is easy to show that the levels of the functions are not even 
always connected. There are examples such that when interpolating between a pair of moving walkways a, b 
and a', b', the travel time starts to increase, then decrease afterwards because the customer uses the walkway 
in the other direction. 

Remark 2. In the above proof, we implicitly use the assumption on the computation model that constant- 
size subproblems can be solved in constant time. This is a standard assumption, but it can make the given 
algorithms difficult to implement in practice. 

4.2 Travel time diameter 

As a next step toward an algorithm for minimum travel time diameter location, we give an 0(nlog n) time 
decision procedure for the travel time diameter given a moving walkway (not necessarily horizontal). Note 



that the travel time diameter for a moving walkway in the plane generalizes the travel time diameter for a 
moving walkway on a line, so the 0(n log n) lower bound given in theorem 2 still holds and our algorithm 
is optimal. 

Problem 4 Given a set P of n points, a real number v > 1, and a pair of points a, b, find 

max dt; me (s, t, a, b). 

s,t£P 

Lemma 7. The decision version of Problem 4 can be solved in optimal 0(n log n) deterministic time. 

Proof. The decision problem is to check, for a given walkway a, b, whether the largest travel time 
max s t6 p <itime(s, a : b) between two points of P is not greater than a given value y. 

We first observe that the pairs of points s,t G P such that d%i me (s, t, a, b) < d(s, t) are all separated 
by the bisector of the line segment [ab]. If both points are on the same side of the bisector, then it is always 
more advantageous to go directly from s to t than to use the walkway. So as a first step, we can partition the 
point set in two subsets R = {s E P : d(s, a) < d(s, b)} and B = {t E P : d(t, a) > d(t, &)}, say of red 
and blue points, according to which side of the bisector the points lie. Then we have : 

max dti mc (s, t, a, b) = maxjmax d(s, t), max d(s, t), max c? timc (s, i, a, 6)}. 

s,tEP s,t£LR s,t£B s£R,t£B 

So by first computing, in time 0(n log n), the Euclidean diameter of the red and blue sets, the travel time 
diameter problem boils down to computing the red-blue travel time diameter. 



y - d(s, a) - 



Fig. 5. The travel time disk B(s, y). 



The red-blue travel time diameter is not greater than y if and only if each point t € B is contained in 
the intersection of the travel time disks of the points in R, where the travel time disk B(s, y) of s 6 R is 
the set of points reachable in time y from s. We observe that for a red-blue pair s S R, t G B, we have 
dtime(s, t, a, b) = mm{d(s, t),d(s, a) + -d(a, b) + d(b, t)}, which means that for going from a red to a blue 
point, the walkway is used from a to b. So the set B(s, y) for a point s E R is either (see figure 5): 

- the disk of radius y centered at s if y < d(s, a) + -d(a, b), or 

- the union of a disk of radius y centered at s and a disk of radius y — d(s, a) — \d{a, b) centered at b, 
otherwise. 



To check the intersection condition, we proceed by first sorting the points in R in nondecreasing order 
of their distance to a. We denote by s\, S2, ■ ■ ■ , S|.r| the points in R, such that d(s\, a) < d(s 2 ,a) < . . . < 
d(s\ R \,a). 

For each point t e B, we identify, in 0(log \ R\) time, the smallest index i such that Vj > i : d(sj,a) > 
y — ^d{a 1 b) — d{b 1 1). All the red points Sj for j > i are such that we cannot go from Sj to t in time less 
than y by using the walkway. So for all these points, we must be able to reach t directly, which means that t 
must lie in the intersection of the Euclidean unit disks T>(sj , y) with centers Sj and radius y. 

To check the latter condition, we can use a simple point location structure for disk intersections. The 
difficulty lies in the fact that we must keep a version of the point location structure for each intersection of 
the form Cj = Plj=i ^( s j) v) f° r * = 1, 2, . . . , To this purpose, we construct the intersections Cj by 
adding the disks one by one in an incremental data structure. This incremental data structure can then be 
made persistent using known techniques [13], with a O(l) multiplicative cost. 

The structure we use is made of two balanced search trees, one for the upper boundary of the intersection 
Cj, and another for the lower boundary. These trees allow to compute, in time logarithmic in the number 
of vertices of Cj, for any value of the x-coordinate, the highest (resp. lowest) point on the boundary of C, 
having this x-coordinate. From this, we can check if a point belongs to Cj. We now need a technical lemma. 

Lemma 8. Consider the intersection I ofi unit disks, and an additional, distinct, unit disk V. The boundary 
ofV intersects the boundary of I in at most two points. 

This shows that the number of vertices of Cj is 0(i) and that the point location data structure answers queries 
in time O(logi) = O(logn). Using the partial-persistence transform [13], the additional space cost due to 
persistence is the total number of pointer changes on the structure over a sequence of operations. Since there 
are at most two new vertices on the intersection of the disks for each new value of i, the total number of 
changes is 0(|i?|). So the persistent data structure uses linear space. 

So for each point t e B, we can check in 0(log n) time whether ted for the computed index i. Hence 
the whole algorithm takes 0(n log n) deterministic time. □ 

4.3 Minimum travel time diameter location of horizontal moving walkways 

Now given our algorithm for computing the travel time diameter and the observation that the location for n 
source-destination pairs is LP-type, we are able to solve the minimum diameter location problem, which is 
in fact an implicit quasiconvex program. 

Problem 5 Given a set P of n points, and a real number v > 1, find a horizontal moving walkway [a* , b*] 
of speed v minimizing the travel time diameter of P: 

min max dtimo(s, t, a, b). 

a,b:a y =b y s,t£P 

Theorem 6. Problem 5 can be solved in 0(n log n) randomized expected time. 

Proof. The result is by using lemma 5. This requires to be able to decompose a set of points P into r groups 
Pi, each of size at most a\P\, such that the set of constraints encoded by P (the pairs s,t € P) is the union 
of the corresponding sets for the Pj. Here, we can partition the set P in three equal-sized subsets, say Q, R 
and S. We form the following three groups of points : Pi = Q U R, P 2 = R U S and P 3 = Q U S, each 
group containing 2n/3 points. We have r = 3 and a = 2/3. This decomposition satisfies requirement 3 of 
the lemma, because each pair of points is included in at least one group. 



Another thing we need from requirement 2 of the lemma is an algorithm for the satisfaction- violation 
test for a group of constraints. This corresponds to a decision algorithm for the travel time diameter, which, 
from the previous lemma, can run in 0(n log n) time. 

Hence from lemma 5, the optimization problem can be solved within the same time bounds, in the 
randomized expected sense. □ 



4.4 An algorithm for moving walkway location in the plane 

So far we only considered the constrained problem in which the moving walkway must be horizontal. The 
unconstrained problem is as follows. 

Problem 6 Given a set P of n points and a real number v > 1, find a moving walkway [a* , b*] of speed v 
minimizing the travel time diameter of P: 

min max dtimcfs, t, a, b). 

a.b s,teP 

It is easy to verify that the constraints are not quasiconvex anymore, so we are not able to apply the 
previous randomized techniques. We can solve the problem approximately, however, using the following 
trick: rotate the set of points in 0(l/e) directions, and for each direction solve the constrained version of the 
problem. 

Theorem 7. A ( 1 +e) -approximate solution for problem 6 can be found inO(^n log n) randomized expected 
time. 

Proof. For each angle ie/v, with i = 0,1,..., [27nj/eJ , we rotate the set of points of this angle and solve 
Problem 5. Among all the solutions computed, the one minimizing the diameter is retained. 

The angle of the walkway returned by this procedure and the angle of an optimal moving walkway a* , b* 
differs by at most a = e/v. Also, the returned solution is not worse than a walkway a, b with the same 
orientation and such that a = a* and d(a, b) = d(a*,b*). Then for any pair s, t we have 

d{s, a) + -d{a, b) + d(b, t) < d{s, a*) + -d(a*,b*) + d(b*,t) + d(b, b*) 

< d (s,an + W, b n + d(b*, t ) + 2 S in a d(a*,b*) 
v 2 

< d(s,a*) + -d(a* ,b*) + d(b* ,t) + ad(a* ,b*) 

v 

< d(s,a*) + -d(a*,b*)+d(b*,t) + £ d(a*,b*) 

v v 

< (l + £)(d(s,a*) + -d(a*,b*) + d(b*,t)). 

v 

Note that the dependence on v is necessary. In particular, if v — > +oo, the above procedure is not applicable. 

□ 

Remark 3. By combining the above trick with a grid technique (see for instance Har-Peled [17]), we can 
obtain a location algorithm whose complexity is of the form 0(n + poly(e^ 1 , v)). The main idea is to snap 
the points in P to a grid of size e -1 x e _1 , and work on the point set formed by taking the centers of all 
nonempty cells. We can then run the above algorithm for n = 0(s~ 2 ). 



Remark 4. Finding an efficient exact algorithm for problem 6 is the natural extension of this work. It is easy 
to see that the problem is polynomial. In the unconstrained setting, every pair of points of the set defines 
a constraint consisting of pieces of 4-variate functions in 5-space. The pieces are algebraic functions of 
bounded degree. Thus, the complexity of the arrangement is polynomial, as well as the time to compute the 
upper envelope (see, e.g., [21]) and to find the minimum value over the envelope. Note however that the 
complexity of the arrangement is in 0(n 8 ), which means that this approach is probably not tight. Whether a 
better upper bound can be derived is left open. 

5 Escalators and Elevators 

In this section, we consider a number of variants of the moving walkway model, by combining a number of 
restrictions on the way the walkway is used. One of the model we obtain is the elevator model, for which the 
location algorithm becomes linear. In the last part of this section, we give a simple algorithm for computing 
the travel time diameter of a set of points with respect to k elevators. 

5.1 Variants 

We consider three restrictions to our model. 

Unidirectional walkways. In that case, we assume that the walkway can only be used for going from a 
to b, and not for going from b to a. If we simply assume this, the diameter of a point set is trivially the 
geometric diameter, since for each pair of points, the walkway cannot be useful in both directions, so the 
travel time in one of the two directions will be equal to the geometric distance between the points. To make 
the model interesting, we must further assume that the travel direction for each pair is prescribed. This can 
be achieved for instance by dividing the point set in two subsets, say red and blue points, and assuming 
that we always travel from a red point to a blue point. In this model, the constraints defining the problem 
are all quasiconvex, whatever the orientation of the walkway is. And since the diameter algorithm can be 
adapted, we immediately have a 0(n log n) randomized expected time algorithm for the optimal location of 
a (unidirectional, red-blue) walkway of arbitrary orientation. 

Mandatory use. With this restriction, we assume that customers are obliged to use the walkway. In that case, 
the algorithm for computing the diameter boils down to computing the two points (one in each color set) that 
are the furthest away from a and b, respectively. The whole location algorithm therefore becomes 0(n) for 
a horizontal walkway. If combined with the previous restriction, we obtain a 0(n) algorithm for the optimal 
location of a (mandatory, unidirectional, red-blue) walkway of arbitrary orientation. We call this model the 
escalator model, as the red and blue subsets can represent source and destination points in two different 
floors of a building. 

Elevators. The elevator model is obtained by assuming the previous two restrictions and letting a = b. In 
this case, the problem is similar to a point location problem mentioned by Eppstein ([14], section 3.2), and 
solvable in 0(n) expected time using Chan's method again. 

5.2 Computing the travel time diameter with respect to fc elevators 

We now give an algorithm that, given the location of k elevators, efficiently computes the diametral source- 
destination pair. This can be considered as a first step towards an algorithm that would locate k elevators 
instead of one. 



Problem 7 Given a set R of red points, a set B of blue points, and a set {ej}j = i ; 2,...,fc of points (elevators), 
find a pair r, b solving the following problem: 

max mm d(r, a) + d(ei,b). 

r£R,b£B i 

We assume that \R\ = \B\ = n. 
Theorem 8. Problem 7 can be solved in 0(kn log*" 1 n) deterministic time. 

Proof. Suppose that R = {ri\i—i t 2,...,n an d B = {&i}i=i,2,...,n- Elevator e\ is used to go from rj to bj iff: 

d(r h ei) + d(e u bj) < d(n, e e ) + d{e e , bj) W = 1, 2, . . . , k (1) 
d{n, ei) - d{n, e e ) < d{e e , bj) - d{e u bj) W = 1, 2, . . . , k. (2) 

V V- ' v ' 

_p(i) 

After having computed the vectors p(*) and Q^) for each and each &j respectively, finding the maximum 
time distance between two points using elevator e\ amounts to finding the maximum distance d(ji, e\) + 
d(ei, bj) between all pairs r^, bj for which dominates pW. This can be done in 0(n log* -1 n) using a 
divide-and-conquer algorithm. Doing this once for each elevator yields the result. □ 
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